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O |. Abstract. The infrared behavior of gluon and ghost propagators in Landau gauge Yang-Mills 

<1} theory has been at the center of an intense debate over the last decade. Different solutions of 

' the Dyson-Schwinger equations show a different behavior of the propagators in the infrared: 

in the so-called scaling solutions both propagators follow a power law, while in the decoupling 
solutions the gluon propagator shows a massive behavior. The latest lattice results favor the 
decoupling solutions. In this contribution, after giving a brief overview of the present status 
of analytical and semi-analytical approaches to the infrared regime of Landau gauge Yang- 
CTN ' Mills theory, we will show how Callan-Symanzik renormalization group equations in an epsilon 

expansion reproduce both types of solutions and single out the decoupling solutions as the 
infrared-stable ones for space-time dimensions greater than two, in agreement with the lattice 
' calculations. 
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1. Introduction 

The infrared regime of quantum chrome-dynamics is one of the prime examples of non- 
perturbative physics giving rise to such fascinating phenomena as the confinement of quarks and 
the dynamical breaking of chiral symmetry. The only method so far to obtain quantitatively 
reliable results directly from the underlying theory has been lattice gauge theory, and it has 
done so only in a limited fashion. For a long time, it was generally believed that analytical tools 
could not be of much help in the quantitative description of the aforementioned phenomena. This 
point of view is now slowly changing. In this contribution, before describing our own approach 
to the subject to which the title refers, we will try to give an overview, from our perspective, of 
what has been achieved so far in (the pure gauge sector of) infrared quantum chromodynamics 
with analytical and semi-analytical methods. Needless to say, progress with these methods has 
also been limited, but there have been some recent very promising developments. 

The story begins in 1978 with the publication of Gribov's seminal work on the quantization 
of non-abelian gauge theories [1]. It was clear by that time that for the quantization of a 
gauge theory in the continuum the gauge needs to be fixed, and that, at least for perturbation 
theory, covariant gauges are particularly convenient. Specializing to the Landau gauge, the 
functional integral over the gauge field then includes two characteristic pieces: a functional 
delta distribution that implements the gauge condition d^Af^ = 0, and a functional determinant 
(Jacobian) that arises when one changes the integration variables from the general gauge fields 
to gauge equivalent fields A^ that fulfill the gauge condition, and the corresponding gauge 



transformations U, A® = (A u )®. The delta distribution can conveniently be written with the 
help of a bosonic auxiliary field called the Nakanishi-Lautrup field, 

S{d lx A%) = j D[B] exp {-i j d D xB a d,A a )j , (1) 

while the Jacobian is represented in a local way as a functional integral over (fermionic scalar) 
ghost and antighost fields, 



det(-^D^ b ) = j D[c,c] exp (J d D x^ 1 d fl Dfc b 



c b ■ (2) 



We are generally working in D-dimensional Euclidean space-time, and the covariant derivative 
(in the adjoint representation) is defined as 

Df = S a %-gf abc Al. (3) 

In the functional integral of the theory (for gauge invariant operators) the integral over the gauge 
transformations U then factorizes and can consistently be omitted, and one is left in the pure 
gauge sector with an integral over A" B a , c a and c a , and the (negative of the) Faddeev-Popov 
action 



S FP = J d D x Qi^i^ + WdpAl + d^D^c^j 



(4) 



in the exponent. 

The original gauge symmetry of the Yang-Mills theory is manifestly broken by the gauge 
fixing procedure, however, the integration measure and the Faddeev-Popov action are invariant 
under the global fermionic BRST transformation s defined as 



fabjb 



sB a = 0, 



sc a = ^gf abc c b c c , sc a = iB a . (5) 

An important property of this transformation (including the transformation of the matter fields) 
is its nilpotency, s 2 = 0. Note that the transformation of the gluon field has the form of an 
infinitesimal gauge transformation with the ghost field replacing the generator of the gauge 
transformation. The Ward identities corresponding to the BRST symmetry are the Slavnov- 
Taylor identities for the correlation functions. 

As it turns out, the beta function for this theory is negative giving rise to the celebrated 
asymptotic freedom of non-abelian gauge theories and the applicability of perturbation theory 
in the ultraviolet regime. In the language of renormalization group theory, non-abelian gauge 
theories have a trivial fixed point (with vanishing renormalized coupling constant <7i?(/u)) which 
is ultraviolet-attractive. On the other hand, the renormalized coupling constant runs into a 
Landau pole at the characteristic scale Aqcd, and perturbation theory (around the trivial fixed 
point) becomes useless in the infrared, i.e., for energy-momentum scales of the order of and 
below Aqcd- 



2. Gauge copies 

In 1978, Gribov [1] observed that the gauge fixing procedure leading to the Faddeev-Popov 
action as described above is actually not consistent in non-abelian theories. The problem arises 
in the change of variables in the functional integral over the gauge field: the choice of a gauge- 
fixed field d^A® = 0, that lies on the same gauge orbit as the initial field A® = (AF) a ^ 



for an appropriate gauge transformation U, is not unique. In other words, there exist different 
gauge fields, called gauge copies or Gribov copies, that fulfill the gauge condition and arc, 
nevertheless, connected by a gauge transformation. In geometric language, the gauge orbits 
intersect the hypersurface defined by d^A^ = more than once. 

In particular, if we look at an infinitesimal gauge transformation generated by oj a , we have 

8^)1 = 8^ + dpD*u/>. (6) 

Then it is possible to have both, A^ and fulfill the gauge condition if oo a is a zero mode of 

the Faddeev-Popov operator (— d^D^). Hence Gribov's proposal to restrict the integration over 
the gauge field to the region where the Faddeev-Popov operator is positive definite, now called 
the (first) Gribov region and denoted by ft. Note that at the boundary of the Gribov region, the 
(first) Gribov horizon dfl, we have that det(— d^D^) = because the Faddeev-Popov operator 
develops a zero mode there. 

Although gauge copies which are related by an infinitesimal gauge transformation are 
excluded by the restriction to ft, it is not clear that finitely related gauge copies cannot exist 
within ft. In fact, it is now known that such copies do exist in ft, and the consequent thing to do 
would be to restrict the integration over the gauge fields further to a smaller region that is called 
the fundamental modular region A. However, the restriction to A is very hard to implement, 
in numerical calculations (lattice gauge theory with Landau gauge fixing) and even more so in 
analytical (or semi- analytical) calculations. Fortunately, according to an argument by Zwanziger 
[2], the difference between the restrictions to ft or to A has no effect on the correlation functions. 

For the rest of this contribution, we will hence consider the restriction of the integral over 
the gauge field to the Gribov region ft. One of the important implications of this restriction as 
compared to an unrestricted integral over the gauge field, according to an argument by Gribov 
[1], is the impossibility of a Landau pole for the renormalized coupling constant to occur (at 
a momentum scale different from zero). Realizing the restriction to ft in an approximate way, 
Gribov was furthermore able to predict the infrared behavior of the gluon and ghost propagators. 
Here we will present an exact way to implement the restriction of the integral over the gauge fields 
which is due to Zwanziger [3] and based on Gribov's ideas, and come back to the propagators 
later. 

Zwanziger found a way to realize the (Heavyside) theta functional for the restriction of the 
integral over the gluon field to the Gribov region as a Gaussian (in the thermodynamic limit): 

J^D[A] = J D[A]9(-d,Df) = J D[A] exp (- 7 4 | d D xh(x)j , (7) 

with the horizon function h(x) defined as 

h(x) = g 2 f abc f cde Al [(-d.D,)- 1 ] ad Al . (8) 
For the correct normalization, the horizon condition 

(h(x)) = 4(N 2 ~ 1) (9) 

has to be fulfilled (in an SU(iV) gauge theory). This condition fixes the new parameter 7 with 
the dimension of mass (in terms of Aqcd ; or vice versa). 

Now h(x) is a complicated nonlocal functional because it contains the inverse of the Faddeev- 
Popov operator. A local form of the theory, necessary for the application of the usual quantum 
field theoretical machinery, can be achieved by introducing two additional auxiliary vector fields 
(and the corresponding antifields), one of them bosonic <p ab and the other fermionic co® b , each 



carrying two color indices (in the adjoint representation). The gauge sector of the theory 
including the restriction to the first Gribov region is then described by a functional integral 
over A®, B a , c a , c a , <#f , (j) ab , L0® b and tiff, and the exponential of the (negative of the) Gribov- 
Zwanziger action [4] 

Sgz = S FP - J d D x (d^Df^ - d^fDfufr + gr hc d^fD b ; C ^ c u d ) 

- 7 2 j d D x (gf abc A^ + gf^Affi + 4(iV 2 - 1) 7 2 ) , (10) 

see eq. (4). The horizon condition can be rewritten in terms of the auxiliary fields as 

{9f abc A^ + ^ c )> = 8(N 2 - 1) 7 2 . (11) 
The Gribov-Zwanziger action leads to the tree-level gluon propagator (in momentum space) 

(A^AK-q)) = -j^j (V - 5 a \2 7 r) D 6(p - q) , (12) 

with the characteristic scale A defined as 

A 4 = 2g 2 N 1 A . (13) 

For p 2 3> A 2 , one recovers the usual (tree- level) ultraviolet behavior, while the propagator is 
heavily suppressed for p 2 <C A 2 and even tends to zero in the limit p 2 — > 0. The tree-level ghost 
propagator, on the other hand, behaves like 

(c a (p)J(-q)) cx 1 5 ab (2ir) D 5(p - q) (14) 

in the infrared regime p 2 <C A 2 . The divergence in the limit p 2 — > is much stronger than for 
the tree-level ghost propagator in Faddeev-Popov theory that does not take the restriction to 
the Gribov region into account. It is clear, then, that the effect of the restriction to the Gribov 
region is rather drastic for the infrared behavior of the propagators. The forms (12) and (14) of 
the propagators have already been found by Gribov in his approximate description [1]. 

3. Dyson- Schwinger equations 

In an initially independent development, the Dyson-Schwinger equations of quantum 
chromodynamics were approximately solved in the infrared regime. As in the previous section, 
we here focus on the pure gauge sector, i.e., Yang- Mills theory. The Dyson-Schwinger equations 
are derived from the Faddeev-Popov action (4), and the restriction to the Gribov region is not 
implemented in any way. While this is clear from the independence of this development, the 
reader is justified to wonder why it would make any sense to present the solution of Dyson- 
Schwinger equations that ignore the restriction to the Gribov region, when we have argued at 
the end of the last section that the restriction to the Gribov region has a severe effect on the 
infrared behavior of the propagators. 

However, as Zwanziger has clarified [5], the restriction to the Gribov region actually leaves 
the Dyson-Schwinger equations unchanged: Dyson-Schwinger equations are derived from the 
fact that the functional integral over a (functional) derivative vanishes. When one restricts 
the functional integral over the gauge field in the Faddeev-Popov formulation to the Gribov 
region Q, contributions from the boundary dQ are to be expected. However, there are no such 
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Figure 1. Simplified Dyson- Schwinger equations for the gluon and ghost propagators containing 
the dominant contributions in the infrared regime. 

contributions in this case since the integrand vanishes on 90 due to the presence of the Faddeev- 
Popov determinant as remarked earlier, and the resulting Dyson-Schwinger equations are the 
same as for the unrestricted functional integral. 

Solutions to the Dyson-Schwinger equations for the gluon and ghost propagators have been 
found, in a certain approximation, first in refs. [6, 7]. The results are an infrared suppressed 
gluon propagator and an infrared enhanced ghost propagator as compared to the tree-level 
propagators of Faddeev-Popov theory, which is qualitatively consistent with the results (12) and 
(14) of the Gribov-Zwanziger framework. Then, with hindsight, we can simplify the Dyson- 
Schwinger equations in the infrared regime by omitting the diagrams that contain gluon loops, 
or in general, in the Dyson-Schwinger equation for a given correlation function only keep the 
diagrams with the biggest number of ghost loops since the latter dominate the infrared behavior. 
The dominance of the diagrams with the biggest number of internal ghost propagators in the 
infrared regime is sometimes referred to as ghost dominance. 

Furthermore, Taylor's non-renormalization theorem [8] allows us to replace the full ghost- 
gluon vertex by the bare one. That this replacement is valid to a good approximation even in the 
non-perturbative regime has been checked in lattice calculations [9]. With these simplifications, 
the Dyson-Schwinger equations for the gluon and ghost propagators take the form depicted in 
fig. 1. This simplified system of equations not only correctly describes the extreme infrared limit 
of the theory, but it also gives a qualitative description of the extreme ultraviolet, asymptotically 
free regime, and it can thus serve as a useful laboratory for the study of several properties of 
the complete theory [10]. For a quantitatively correct description of the ultraviolet regime, in 
particular the correct value of the (one-loop) beta function, however, one needs to take additional 
contributions into account in the first equation of fig. 1, and more sophisticated approximation 
schemes have been devised to this end [11, 12]. 

If one is exclusively interested in the extreme infrared regime, the simplified Dyson-Schwinger 
equations can even be solved analytically. In order to describe the solutions, we will introduce 
some notation: we write the propagators in terms of their dressing functions G(p 2 ) and F(p 2 ), 



For the analytical solutions, we make power-like ansaetze for the dressing functions (this behavior 




(c a (p)c b (-q)) 



HP 2 ) xab 



5 ab (27r) D 5{p-q). 



(15) 



is familiar from the theory of critical phenomena): 

G{p 2 ) oc {p 2 y aG , F{p 2 ) oc (p 2 )- aF . (16) 

Inserting the ansaetze (16) in the system of approximate Dyson- Schwinger equations in fig. 
1, one obtains, for any dimension D between two and four, two consistent solutions known as 
scaling solutions [5, 13]. One of these solutions is 

a F (D) = 1 ^, aG{D ) = - I l (17) 

in dimension D, for 2 < D < 4 (this solution is not entirely consistent in dimensions two and four 
due to logarithmic corrections). The infrared exponents of the other scaling solution have to be 
obtained, for general D, from the solution of a transcendental equation. A good approximation 
(within 2% error for any D with 2 < D < 4) is given, however, by the simple formulae 

, , D-l . , 16 -D 
«fP) = — > a G (D) = — . (18) 

For both solutions (also the exact version of the second solution), the infrared exponents 
fulfill the sum rule 

D-4 

a G + 2a F = — ^— (19) 

in any dimension D. The solutions also have in common that ap > and 1 + olq < (for 
D > 2 in the first solution), implying that the ghost propagator diverges more strongly in the 
infrared than the (Faddeev-Popov) tree-level propagator, while the gluon propagator is strongly 
suppressed in the infrared, i.e., it tends to zero in the limit p 2 — > 0, just like the tree-level 
propagators in the Gribov-Zwanziger framework. 

Taking into account that the Dyson-Schwinger equations are not altered by the restriction 
of the functional integral to the Gribov region as argued at the beginning of this section, we 
interpret the solutions (17) and (18) as (approximately) incorporating the quantum corrections 
to the tree-level results (12) and (14) presented in the previous section. It is also clear that 
the power-like behavior (16) of the solutions can only be obtained by summing up the quantum 
corrections (in some approximation) to any loop order. On the other hand, it is a priori not so 
clear which of the two solutions (17) and (18) is the physically relevant one for general dimension 
D. Also, one would like to have a systematic scheme to improve on the approximations employed 
in the derivation of the Dyson-Schwinger equations in fig. 1, and it has turned out to be very 
difficult to formulate such a scheme. 

An important consequence of Taylor's non-renormalization theorem is the natural definition 
of a running coupling constant as 

9Up 2 ) = (p 2 ) {D ~ m G(p 2 )F 2 (p 2 )g 2 . (20) 

The factor (p 2 )( D ~ 4 )/ 2 compensates for the dimension of the bare coupling constant g 2 (the 
dressing functions G(p 2 ) and F(p 2 ) are by definition dimensionless), so that g\{p 2 ) is a 
dimensionless quantity. Then it is clear from the sum rule (19) that the running coupling 
constant tends to a finite fixed point value in the infrared limit. This value can be determined 
analytically for the two scaling solutions (17) and (18), and it turns out to vary continuously 
with the dimension D. As an example, we give the (exact) fixed point value in D = 4 for the 
second solution (18), 

^ = 8.915..., (21) 

and mention that g\ tends to zero in the limit D — > 2 for the first solution (17). 

Finally, we mention that the same results have also been obtained (in D = 4 dimensions) with 
a different technique, the functional renormalization group equations [12, 14] which, however, 
have many features in common with the Dyson-Schwinger equations. 



4. Decoupling solutions 

One of the characteristic features of the scaling solutions (including their tree-level forms (12), 
(14)) is the infrared enhancement of the ghost propagator. In fact, it may be shown [4] that the 
horizon condition (9) is equivalent to 

lim F^ 1 (p 2 ) = 0, (22) 

the infrared divergence of the ghost dressing function. It is well understood by now that 
implementing the condition (22) necessarily leads to the scaling solutions, be it in the tree- 
level form (12) and (14) or in the "resummed" form (17) or (18). 

There is another important aspect of the Gribov-Zwanziger framework [15]: the horizon 
condition is identical to the Kugo-Ojima confinement criterion [16, 17] which guarantees (in 
Faddeev-Popov theory) the existence of a global color charge that is BRST-exact and hence 
implies confinement in the sense that all physical states are color singlets. Kugo and Ojima 
construct the physical Hilbert space of non-abelian gauge theories starting from BRST invariance 
and show the unitarity of the S'-matrix [16]. The equivalence of the confinement criterion (and 
hence the horizon condition) to the infrared divergence (22) of the ghost dressing function in 
the Landau gauge has been shown in this framework, too [17]. 

At first sight, everything seems to fit together nicely: the restriction of the functional integral 
over the gluon field to the Gribov region, necessary because of the existence of gauge copies in 
non-abelian gauge theories, leads to the horizon condition or the Kugo-Ojima criterion and hence 
to the confinement of color charges. However, there is a fundamental flaw in this argument: the 
BRST symmetry of the Faddeev-Popov action can be extended to the Gribov-Zwanziger action 
(10) only for 7 2 = 0, in which case the integration over the gluon field is effectively unrestricted, 
see eq. (7). In other words, the restriction to the Gribov region breaks the BRST invariance on 
which the whole Kugo-Ojima construction rests. 

There is a simple intuitive argument that shows how the breakdown of BRST invariance is 
an inevitable consequence of the restriction to the Gribov region [18]: the BRST transformation 
of the gluon field has the form of an infinitesimal gauge transformation, as already remarked 
following eq. (5). But close to the Gribov horizon, an infinitesimal gauge transformation can 
relate gauge copies, see eq. (6), and hence the transformed field falls outside the Gribov region. 

The violation of BRST symmetry is not in itself a catastrophe for the theory. Many of us 
have grown up thinking of gauge invariance and BRST invariance as synonyms, mostly because 
the (global) BRST symmetry implies the Slavnov- Taylor identities. However, this is really only 
true for the Faddeev-Popov action and an unrestricted functional integral over the gluon field. If 
we restrict the intregration to the Gribov region, as it turns out to be necessary in order to avoid 
gauge copies, the (gauge fixed) quantum theory is not invariant under BRST transformations 
any more. Of course, this is a serious problem for the construction of the Hilbert space of 
the theory, just because the standard construction (of Kugo and Ojima) is based on BRST 
symmetry. A solution to this problem has not yet been devised. 

But there is yet another twist to our story. With the analytical and semi-analytical predictions 
of the propagators in the infrared regime at hand, it seemed natural to corroborate these 
predictions with the help of completely non-perturbative lattice calculations. While it is not 
necessary (nor natural) to fix the gauge in these calculations, it is also not too difficult to 
implement the Landau gauge on the lattice: for a given gauge field configuration A a ^ consider 
the functional 

J d D x(A u )°(A u )° (23) 

for all possible gauge transformations U, i.e., along the gauge orbit that contains A®. Then the 
extrema (or saddle points) of this functional along the orbit correspond to the fields A a = {A u ) a 



that fulfill the gauge condition d^A^ = 0, the local minima correspond to fields A® in the first 
Gribov region Q, and the absolute minima to fields in the fundamental modular region A. For 
the huge lattices that are necessary to measure the infrared behavior of the propagators, it is 
extremely difficult (numerically costly) to determine the absolute minimum of the functional 
(23) along a given gauge orbit, however, it is relatively easy to find a local minimum and hence 
to restrict the gauge fields to the Gribov region. 

The first lattice results for the propagators in the Landau gauge (with the restriction to 
the Gribov region U) show a clear increase of the ghost dressing function F(p 2 ) towards the 
infrared, and a maximum and subsequent marked decrease of the gluon dressing function G{p 2 ) 
as one follows p 2 from the ultraviolet to the infrared, in agreement with the scaling solutions 
of the Dyson-Schwinger equations. However, it was difficult to discern for a long time whether 
the ghost dressing function would actually diverge in the limit p 2 — >■ and whether the gluon 
propagator itself would vanish in the same limit. 

Meanwhile, it was discovered that the Dyson-Schwinger equations of the theory permit 
another type of solutions now known as decoupling solutions [19, 20, 21] (see also [12]). In these 
solutions, the ghost propagator is infrared enhanced only by a finite constant, i.e., the ghost 
dressing functions tends toward a finite constant in the limit p 2 — >■ 0, while the gluon propagator 
has a finite (nonzero) limiting value similarly to the propagator of a massive particle. Then 
one still has ghost dominance in the sense that the gluon propagator is infrared suppressed in 
comparison with the ghost propagator, and the reduction of the full Dyson-Schwinger equations 
to the approximation represented in fig. 1 is justified in the infrared regime. 

In terms of the infrared exponents defined in eq. (16), one has for the decoupling solutions 
that 

«F = , ac = -1 (24) 

in any dimension D (coincident with the first scaling solution (17) in the limit D — > 2). In 
particular, the decoupling solutions do not respect the sum rule (19) for the exponents (except 
at D = 2). Furthermore, the running coupling constant defined in eq. (20) tends to zero in the 
limit p 2 — > for the decoupling solutions (for dimensions D > 2). 

The discovery of the decoupling solutions of the Dyson-Schwinger equations initiated a long 
and often passionate debate in the community about the true infrared behavior of the Landau 
gauge propagators. In 2007, at the lattice conference in Regensburg, three lattice groups 
presented their results on the infrared behavior of the propagators in Landau gauge Yang- 
Mills theory on very large lattices [22, 23, 24]. The majority of the community has interpreted 
these results as confirming the decoupling solutions of the Dyson-Schwinger equations in three 
and four space-time dimensions. In the same year, lattice calculations in two dimensions [25] 
confirmed the realization of the second scaling solution (18). 

At first sight, the decoupling solutions are at odds with the tree-level propagators (12) 
and (14) and the horizon condition in the form (22) that result in the Gribov-Zwanziger 
scenario. However, indications were found in ref. [18] (see also [26]) that the auxiliary fields 
form a condensate (</>jf <^ 6 ) (among other possible condensates). As a consequence, the horizon 
condition (11) does not imply the infrared divergence (22) of the ghost dressing function any 
more, and through the linear coupling of the fields and <#f to the gluon field, see eq. (10), 
the condensate leads to a mass term for the gluons. 

In the simplest case [18] in this so-called refined Gribov-Zwanziger framework, i.e., with the 
condensates taken into account, the tree-level propagators are of the form 

(Al(p)AU- q )) = p4 f + ^ 2 2 + A4 (V - Mil) 6« b (2n) D 6(p - g) (25) 



(see eq. (13) for the definition of the characteristic scale A), and 

(c a (p)c*(-q)) cx 1 5 ab {2v) D 5{p - q) (26) 

in the infrared regime. The gluon propagator can also be written in terms of an effective 
momentum-dependent mass M e g(p 2 ), 

p 2 + M 2 _ 1 

p 4 + M 2 p 2 + A 4 - p 2 + M cS (p 2 ) ' 

where 

^b 2 ) = • (28) 

In four space-time dimensions, a successful fit of the tree- level gluon propagator in the refined 
Gribov-Zwanziger framework with an additional (yl^^)-condensate, 

W«4(-«)> = ^p^hmv - ^) a °' (2 " )Da(p - "» ' (29 » 

to the lattice data in the infrared and intermediate momentum regimes has been performed in 
ref. [27] (see also [28]). It is important to clarify that the different mass parameters that have 
appeared so far are not additional parameters of the theory, but can rather be determined, at 
least in principle, as multiples of Aqcd from a calculation of the effective potential [18]. 

At this point, from the perspective of the Dyson-Schwinger equations, the decoupling solution 
as well as both of the scaling solutions are consistent and valid solutions, and from these equations 
alone there is no reason to prefer one of the solutions over the others. In the refined Gribov- 
Zwanziger framework, on the other hand, the decoupling solution with the tree- level propagators 
(25) or (29) and (26) is preferred according to the lowest-order calculation of the effective 
potential for various local composite operators [26]. The problem with this latter approach 
is rather of a technical nature: it appears extremely cumbersome to carry the calculations to 
higher orders. 

We finally mention an important property of the scaling as well as the decoupling solutions: 
all the solutions found so far violate positivity, i.e., the spectral function in the Kallen-Lehmann 
representation of the gluon propagator is not positive semidefinite implying that the gluon cannot 
be a physical (massive or massless) particle (see, e.g., refs. [7, 12, 18, 29, 30]). 

5. Callan-Symanzik equations 

5.1. A gluon mass term 

For the rest of this contribution, we will describe a recent approach to the subject from 
a renormalization group perspective [31]. We start from the Dyson-Schwinger equations of 
Faddeev-Popov theory which, as explained before, do not change when one restricts the 
functional integral over the gluon field to the Gribov region fi. Since the iterative solution 
of the Dyson-Schwinger equations generates the perturbative expansion, the same is true for 
the perturbative series. On the other hand, the (re) normalization conditions that define the 
quantum theory can certainly change as a result of the restriction. 

One of the most important consequences of the restriction to the Gribov region, as mentioned 
in the last section, is the breaking of BRST invariance. The resulting violation of the Slavnov- 
Taylor identites leads, almost inevitably, to the generation of a gluon mass term as the only 
possible properly relevant term in the (effective) action with the field content of Faddeev-Popov 
theory, i.e., without the additional auxiliary fields of the Gribov-Zwanziger action. Note that 



(27) 



a ghost mass term cannot be induced for the same kinematical reasons that lead to Taylor's 
non-renormalization theorem. We will, hence, include a gluon mass term 



Jd D x±A«m 2 A« (30) 



in the Faddeev-Popov action. 

Such an action has been considered first, as a special case among a much more general class 
of actions, in ref. [32] where the extension of the BRST transformation and the Slavnov- Taylor 
identities to the massive case is discussed (however, the Nakanishi-Lautrup field is not introduced 
in ref. [32]). The Faddeev-Popov action with a gluon mass term has also been used in a recent 
perturbative calculation [33, 34] to which we come back later. Since the mass term originates 
from the breaking of the BRST symmetry as a consequence of the restriction to the Gribov 
region, we expect the (renormalized) mass parameter to be of the order of the characteristic 
scale Aqcd- It is then clear that the introduction of the mass term has no important effect in 
the ultraviolet regime, i.e., for momenta p 2 3> Aq CD . In particular, the value of the ultraviolet 
beta function is unchanged and, of course, the theory is asymptotically free. This is consistent 
with the fact that the restriction to the Gribov region is irrelevant for the perturbative expansion 
in the ultraviolet regime where only "small" fields A® contribute. 

Now, we would at least naively expect that in the extreme infrared regime p 2 <C m 2 the mass 
term (30) dominates over the other term in the action that is quadratic in A®, 

J (IrF A X~ P) ( ^ V " P ^ K[P) ■ (31) 

In fact, this is true also for the loop corrections as we will now illustrate with the help of a 
relevant example: consider the one-loop ghost self-energy 
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It is worth mentioning that the tree-level gluon propagator in Landau gauge is exactly transverse 
in momentum space even though the mass term (30) has a longitudinal part. The diagram (32) 
can be evaluated for arbitrary dimension D, e.g., by introducing Feynman parameters. We 
will here concentrate on D = 2, which will turn out to play a particularly important role, and 
comment on other dimensions later. Then, evaluating (32) in D = 2 and expanding the result 
in powers of (p 2 /m 2 ) for the extreme infrared behavior gives 
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On the other hand, neglecting the term (31) in the action corresponds to the replacement 
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for the tree-level gluon propagator. Making this replacement in the loop diagram (32) leads to 
the result 
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(35) 



(plus terms of order e) in D = 2+e dimensions, with an arbitrary unit of mass k. The appearance 
of an ultraviolet divergency in two dimensions is due to the radical change in the ultraviolet 
behavior of the gluon propagator. To be precise, the symbol g in eq. (35) stands for the bare 
coupling constant multiplied with a power of k such that g has the dimension of mass, and it 
coincides with the coupling constant appearing in eq. (33) only in the limit D — > 2. 

The results (33) and (35) are obviously different. However, if we implement an appropriate 
normalization condition for the term proportional to p 2 in the proper ghost two-point function, 
then after renormalization the two results will only differ by contributions of the order (p 2 /m 2 ) 
relative to the leading term, since renormalization takes care of the constants multiplying p 2 . 
Expressions (33) and (35) are then equivalent in the infrared limit, and in this sense the 
replacement (34) of the tree-level gluon propagator is justified. The argument can easily be 
extended to any dimension D in the range 2 < D < 4. 

We can, hence, neglect the term (31) in the action as long as we are only interested in 
describing the extreme infrared regime of the theory. These modifications, the introduction 
of a mass term for the gluons and the subsequent omission of the term (31), have profound 
consequences for the renormalization group anaylisis of the theory. We begin this analysis by 
considering the free (g = 0) part of the action which is now given by 

Sir = / d D x Q Al m 2 A^ + iB a d^ + d^d^ . (36) 

Performing a Wilsonian renormalization group transformation on this free field theory (see, 
e.g., ref. [35]), the only nontrivial step is the rescaling of the space-time coordinates x and the 
fields. Writing the rescaling of the coordinates as x — >■ x/s with a parameter s > 1, the action 
(36) is invariant under the renormalization group transformation if we rescale the gluon field as 

A a {x) ^ s D/2 A a (sx y (37) 

This behavior under rescaling defines the canonical dimension of the gluon field which is hence 
D/2 and manifestly differs from its usual canonical dimension (D — 2) /2. The reason is that the 
action (36) corresponds to the high-temperature fixed point rather than the usually considered 
critical Gaussian fixed point [35], as far as the gluon field is concerned. The high-temperature 
fixed point is naturally related to the model of the stochastic vacuum. The canonical dimension 
of the ghost fields, on the other hand, is the ordinary (D — 2)/2. 

Reintroducing the coupling terms of Faddeev-Popov theory into the action, the scaling 
dimensions around the infrared fixed point action (36) turn out to be different for the different 
couplings. Concretely, both the three-gluon and the four-gluon couplings become perturbatively 
irrelevant, while the ghost-gluon coupling is relevant only for dimensions D < 2 (and marginal 
at D = 2). In other words, the upper critical dimension of the theory is two. This can also be 
concluded from an infrared power counting analysis to all orders of perturbation theory [34]. 

Hence, for an infrared analysis of the theory in dimensions D > 2, we would be justified to 
work with the free action (36) alone, which obviously directly leads to the decoupling solution for 
the propagators (in the extreme infrared limit). It will turn out to be very interesting, however, 
to include the coupling that is marginal at D = 2 dimensions and to treat dimensions above two 
in an epsilon expansion. We will hence add the coupling term 

Sfc = m k^^I 2 gf abc J d D x fygMj c c (38) 

to the free action (36), where we have introduced the dimensionless coupling constant g. 
Incidentally, keeping the ghost-gluon coupling while neglecting the three- and four-gluon 
couplings is precisely equivalent to keeping only the diagrams with the biggest number of 
(internal) ghost propagators, that is, ghost dominance. 



5. 2. Infrared fixed point 

In the theory defined by S® R + Sj R , the one- loop expressions for the propagators are precisely 
those represented in fig. 1, except that the propagators in the loop diagrams should be replaced 
with the tree-level propagators. We have given the result for the ghost self-energy in D = 2 + e 
dimensions already in eq. (35). For the gluon self-energy, we get in D = 2 + e dimensions 
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It is important for the consistency of the renormalization procedure that the counterterm can 
be chosen local, in this case proportional to m 2 (5 M „. The gluon self-energy is not transverse in 
momentum, which is not a surprise since BRST invariance is broken. Only the transverse part, 
given by the coefficient of S^ u in the expression (39), enters the result for the gluon propagator. 

1/2 

For the renormalization of the theory, we introduce renormalized fields as = Z ^ A R 

1/2 

and c a = Z c c R and normalize the renormalized propagators to 
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(40) 



at the renormalization scale fx. Using the expressions (35) and (39) for the self-energies, the 
normalization conditions (40) define the field renormalization constants Za and Z c as functions 
of the scale [i. The one-loop anomalous dimensions are then easily calculated to be 
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We define the dimensionless renormalized coupling constant as usual from the renormalized 
vertex function at the symmetric point, 
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(cf. eq. (38)), where we have omitted the Lorentz and color indices on the symbol T 5 Ac of the 
vertex function. We can now calculate the beta function. Taylor's non-renormalization theorem 
implies that the bare vertex function at the symmetric point is independent of the scale n, so 
that the one-loop beta function becomes 
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From eq. (43) we identify two fixed points for e > (D > 2): a trivial infrared-stable one, 
and an infrared-unstable one at 
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We begin with the discussion of the infrared-unstable fixed point. The ^-independence of the 
bare propagators leads to differential equations for the renormalized propagators: using the 
definition of the anomalous dimensions we obtain the Callan-Symanzik equations 

/i 2 ^ (A% p {p)A b R ^- q )) = - lA (A\ p (p)A b R ^-q)) , 

^ <4(P)4(-S)> = "7c <c&(p)4(-?)> • (45) 

Inserting the fixed point value (44) of the coupling constant in the anomalous dimensions, the 
differential equations for the renormalized propagators may be integrated with the normalization 
conditions (40) as initial conditions. The results are 

„2 \ e/2 



(c R (p)c b R (-q)) = j 2 (^y /2 5 ab (2*) D 8{p - q) . (46) 



This is exactly the infrared behavior of the first scaling solution (17). So the renormalization 
group is able to reproduce one of the scaling solutions in the infrared regime and furthermore 
demonstrates that it corresponds to an infrared- unstable fixed point. 

We should clarify that we refer to eqs. (45) as (simple examples of) Callan-Symanzik equations 
for lack of a better name and following the (rather recent) general use of the term, see, 
e.g., ref. [35]. However, these are not precisely the equations originally derived by Callan 
and Symanzik, but rather homogeneous renormalization group equations that result from the 
reparameterization invariance of the bare correlation functions. 

Now we turn to the trivial infrared-stable fixed point. The fact that the stable fixed point is at 
g R = implies that perturbation theory with the action S® R + Sj R should give a good description 
of the infrared regime. In fact, it was shown in refs. [33, 34] that one-loop perturbation theory 
reproduces the lattice results very well in four dimensions and at least qualitatively in three 
dimensions. The one-loop contributions that include the three- and four-gluon vertices were 
taken into account in this comparison, however, they are suppressed by powers of (p 2 /m 2 ) for 
small momenta. 

Actually, the renormalization group approach can be used to go beyond the mere fact that 
the decoupling solution in the extreme infrared limit corresponds to a stable fixed point. The 
way the coupling constant approaches the fixed point value (zero) can be obtained by integrating 
the differential equation (43) for g~R(n). The result is 
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where we have introduced a new characteristic scale A. It is defined as 
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and should eventually be related to Aqcd (this relation, however, is beyond the scope of the 
present approach). Substituting the result (47) for g in the anomalous dimensions (41), the 



integration of the differential equations (45) yields 



<4(p)4(-«)> = J, ^(^5(p - ,) . (49) 

The infrared behavior of the propagators predicted by eq. (49) is at least qualitatively 
confirmed by lattice calculations in D = 3 dimensions [36]. In D = 4 space-time dimensions 
(e = 2), the increase of the gluon propagator proportional to p 2 starting from a non-zero constant 
at p 2 = as described by eq. (49) is actually seen [37] in lattice calculations at lattice parameter 
(3 = 0, which in continuum language means that the proper Yang-Mills action 

jd D x- A F^ v (50) 

is omitted in the probability distribution, and only the delta functional (1) and the Faddeev- 
Popov determinant (2) (and the restriction to the Gribov region) are left. 

In order to reproduce the D = 4 lattice results with the Yang-Mills action taken into account, 
we have to reintroduce the term (31) into the action S® R + 5j R . The reason is that this term, 
although perturbatively irrelevant given the scaling behavior (37) of the gluon field, is of the 
same order in (p 2 /m 2 ) as the correction to the tree-level gluon propagator obtained from the 
renormalization group improvement in eq. (49) for e = 2. Since at the one-loop level there are no 
quantum corrections to (31) considered as a composite operator, we can just add this term to the 
renormalization-group improved result (49) for the proper gluonic two-point function (inverse 
propagator for the transverse part), which results in 

<^>)A*(-*)> = (p 2 + ™ 2 \t9$) 1 i spa " P ~f) 5ab{2n)D5{p - q) • (51) 

This is precisely of the form (25) of the infrared propagator found in the refined Gribov-Zwanziger 
framework (without an additional (A^ylp-condensate), with the effective momentum-dependent 
mass given by (cf. eqs. (27), (28)) 

M^H^^. (52) 

A different renormalization scheme based on the full Faddeev-Popov action (4) plus the gluon 
mass term (30) has been proposed in ref. [34]. The corresponding Callan-Symanzik equations 
reproduce the decoupling solutions as measured on the lattice reasonably well over the whole 
momentum range in four dimensions, and at least qualitatively in three dimensions. 

In the special case of two space-time dimensions (e = 0), the two fixed points considered so far 
coalesce, and the result is one infrared- unstable trivial fixed point. The instability implies that 
the free action (36) cannot describe the infrared physics in two dimensions. Indeed, calculations 
on the lattice confirm that in two dimensions the decoupling solution corresponding to the action 
(36) is not realized [25]. Instead, apparently the second scaling solution (18) is found, and we 
have not yet succeeded in giving a description of the corresponding fixed point (see, however, 
below) . 



5.3. Lifshitz point 

Up to now, we have not made any attempt to implement the horizon condition in the form (22). 
It is not difficult to do so: on the classical level, for a local (and covariant) formulation, we have 
to replace the ghost term in the action (36) with 

J d D x^d^(-d 2 )d fl c a , (53) 

where we have introduced an arbitrary dimensionful parameter b. The form (53) of the action 
is known in statistical physics as an isotropic Lifshitz point [38, 39]. 

Using the term (53) in the free action changes the canonical dimension of the ghost fields. 
Scale invariance of the free action now implies the scaling behavior 

c a {x) -> S ( D ~ 4 )/ 2 c a (sx) , c a (x) -> S ( D " 4 )/ 2 c a {sx) . (54) 

Then, around the Lifshitz point, the three- and four-gluon couplings are perturbatively irrelevant 
just as before, while the ghost-gluon coupling becomes relevant for dimensions D < 6. 

In order to describe the infrared behavior of the theory with the horizon condition (22) 
implemented in dimensions between two and four, we use an epsilon expansion around D = 6. 
We introduce the coupling term 

S\ K = ™ k^/ 2 gf abc j d D x d^A\ c c (55) 

with the dimensionless coupling constant g (different from the coupling constant in eq. (38)), and 
follow the procedure applied before in close analogy, starting with the calculation of the ghost 
and gluon self-energies, now at dimension D = 6 — e and with the tree-level ghost propagator 
proportional to b 2 /k 4 . We stress that a fine tuning of the tree-level ghost term that appears in 
eq. (36) is necessary in order to fulfill the horizon condition (22) to one-loop order, just like in 
ordinary </> 4 -theory the tree-level mass term has to be fine-tuned to reach the critical point. 

We use normalization conditions and a definition of the renormalized coupling constant ()r{h) 
in strict analogy to eqs. (40) and (42), adapted to the new form (53) of the ghost term and the 
coupling (55). Taylor's non-renormalization theorem continues to hold in the present situation, 
and the result for the one-loop beta function is 

n, - s 2 d 1 _ fe 1 Ng 2 \ , 

^~ 9R) = ^~ 9R = -2~ 9R {2-2Wf) ■ (56) 

It thus turns out that the infrared-stable fixed point for dimensions D < 6 (e > 0) is the 
nontrivial one. The final result for the renormalization-group improved infrared propagators 
corresponding to the infrared-stable fixed point is 
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In the notation of eq. (16), we obtain the following values for the infrared exponents: 



(58) 



The exponents fulfill the sum rule (19). The values (58) should be compared to the (approximate) 
exponents of the second scaling solution (18), 
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The two results coincide (exactly) at D = 6, and in the range 2 < D < 4 they differ the most 
at D = 2 (where the approximation (59) turns out to be exact): ol f = 1/6 vs. a F = 1/5 and 
ac = —8/6 vs. ac = —7/5. However, compared to our experience with the epsilon expansion in 
statistical physics, the coincidence of the infrared exponents (58) and (59) is unexpectedly good, 
and we consider them to be approximations to the same solution of the exact Dyson-Schwinger 
equations of the theory. 

5.4- Discussion 

We have reproduced the two scaling solutions and the decoupling solution of the Dyson- 
Schwinger equations with the help of an analytical renormalization group analysis. Only the 
second scaling solution (18) and the decoupling solution correspond to infrared-stable fixed 
points, while the first scaling solution (17) corresponds to an infrared-unstable fixed point in 
dimensions D between two and four. 

In fact, the Lifshitz point which leads to the second scaling solution is only stable if one 
implements the horizon condition (22) (in addition to the gluon mass term which is a natural 
result of the breaking of BRST invariance). As remarked after eq. (55), the implementation of 
the horizon condition requires a fine tuning of the tree-level term 



a relevant coupling under the scaling behavior (54) of the ghost fields. Consequently, the 
nontrivial fixed point in eq. (56) is infrared-tinstafr/e unless one insists in implementing the 
horizon condition. 

We hence conclude that in dimensions D > 2 the physically realized solution according to 
lattice simulations and the refined Gribov-Zwanziger approach is the only one that corresponds 
to an infrared-stable fixed point. In other words, the only input in the Callan-Symanzik 
approach is the gluon mass term that is generated by the breaking of BRST invariance, 
then the renormalization group flow automatically leads to the decoupling solution without 
further assumptions. In two space-time dimensions, on the other hand, the trivial fixed point 
corresponding to the decoupling solution becomes unstable, and the lattice results point to 
the second scaling solution in this case. However, according to the explanation above the 
fixed point corresponding to this scaling solution should also be unstable. Indeed, calculating 
the one-loop corrections to (60) considered as a composite operator, the fixed point becomes 
even more strongly infrared-repulsive with decreasing dimension. Hence, unless there is some 
reason to implement the horizon condition precisely in two dimensions, we do not yet have a 
satisfactory understanding of the infrared behavior of the theory in two dimensions within the 
Callan-Symanzik approach. 

In conclusion, we have presented a new approach to the infrared regime of Yang-Mills theory in 
Landau gauge that is based on a renormalization group analysis with Callan-Symanzik equations 
in an epsilon expansion. In contrast to the Dyson-Schwinger equations, the Callan-Symanzik 
equations are analytical and inherently systematic as an approximation scheme. Furthermore, 
they are technically relatively simple and economic as compared to the Gribov-Zwanziger 
approach with the additional auxiliary fields necessary for a local formulation. This should 
be an important advantage in the further elaboration of the theory, in particular the inclusion 
of condensates (as in the refined Gribov-Zwanziger framework) and dynamical quarks. 




(60) 
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